Multi-component Transparent Conducting Oxides: Progress in Materials Modelling 
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Transparent conducting oxides (TCOs) play an essential role in modern optoelectronic devices 
through their combination of electrical conductivity and optical transparency. We review recent 
progress in our understanding of multi-component TCOs formed from solid-solutions of ZnO, In2 03, 
Ga2 03 and AI2O3, with a particular emphasis on the contributions of materials modelling, primarily 
based on Density Functional Theory. In particular, we highlight three major results from our 
work: (i) the fundamental principles governing the crystal structures of multi-component oxide 
structures including (Iri203)(ZnO) n , named IZO, and (In203) m (Ga203);(ZnO) n , named IGZO; 
(ii) the relationship between elemental composition and optical and electrical behaviour, including 
valence band alignments; (iii) the high-performance of amorphous oxide semiconductors. From these 
advances, the challenge of the rational design of novel electroceramic materials is discussed. 

PACS numbers: 71.20.-b, 78.20.-e, 71.15.Mb 



INTRODUCTION 

Transparent conducting oxides (TCOs) are denned by 
high electrical conductivity approaching that of a metal- 
lic compound, with the high transmission of photons 
in the visible range of the electromagnetic spectrum. 
The unique combination of these two features makes 
TCOs essential components of modern optoelectronic 
devices [U |2]. The common electronic characteristic of 
n-type TCO materials is a spatially delocalised, low en- 
ergy, and low effective mass conduction band determined 
primarily by the metal s orbitals. Such a state is achieved 
with oxides formed from closed-shell ns° cations such as 
Zn(II), Ga(III) and In(III). 

The prototype n-type TCO is Sn-doped ln 2 3 (ITO), 
which has been the subject of intense study both 
experimentallv(2TfTTj] and theoretically |1 ltfi"5] . Other ex- 
amples of electron conducting TCOs include Al-doped 
ZnO, Ga-doped ZnO and F-doped Sn02- Typically these 
materials have electron carrier concentrations ranging 
from 10 16 - 10 21 cm" 3 , and electron mobility ranging 
from 10 to 1000 cm 2 V _1 s _1 , depending on the material 
quality and growth conditions [2]. 

In contrast to the delocalised conduction band wave- 
functions, the valence band states of hetero-polar metal 
oxides are typically localised O 2p states. Instead of re- 
sulting in mobile electron holes, acceptor doping of these 
materials results in deep localised states, which are sta- 
bilised by local lattice distortions - the formation of small 
polarons[ini [20] • To obtain hole conducting (p-type) 
TCO materials, the principal approach has been to in- 
clude a metal with filled low binding energy states at 
the top of the valence band, such as Cu(I), which can 
facilitate hole transport (21] [22] . Following the work of 
Hosono et al. on CuA10 2 [23 , many ternary Cu oxides 



have been explored for this purpose. The limitation of 
this class of material is the poor hole mobility (< 10 
cm 2 V _1 s _1 ) that arises from the relatively localised Cu 
d bands, as well as the indirect electronic band gaps aris- 
ing from the underlying lattice symmetry. Recent work 
has focused on other delafossite (CuM /// 02, where M 111 
is a trivalent cation) materials including CuGa09|21j. 
CuCr0 2 [21r[2g, CuSc0 2 [27] and CuB0 2 [2H], in addition 
to studies on related ternary copper oxides containing Sr 
and Pb[H[2Hll!2- Unfortunately for all Cu-based mate- 
rials, the poor carrier mobility, which is directly related 
to the parent compound Cu 2 0, is difficult to overcome. 
In some cases, the existence of deep defect centers[3"T]. 
which act as trapping centres for free carriers, pin the 
Fermi level well above the valence band, thus limiting 
the carrier densities achievable by doping. 

The possibility of discovering alternative n-type TCO 
compositions with superior material properties has led 
to the investigation of ternary and quaternary systems, 
with In 2 03(ZnO)„ (IZO) being one particular focus of 
attention [3"2Tf3"7| . These studies have highlighted the im- 
proved chemical and thermal stability of IZO compared 
to ITO, making it more desirable for commercial applica- 
tion. One issue with the utilisation of multi-component 
oxide materials is the difficulty in synthesising high- 
quality crystalline samples at low cost. To overcome this 
limitation, there has been an active interest in growing 
amorphous TCOs (a-TCOs) [35l EHHH] . Amorphous ox- 
ides can exhibit electrical properties comparable to their 
crystalline phases, in contrast to traditional amorphous 
semiconductors. In amorphous systems, the underlying 
crystalline nature and translational symmetry of a mate- 
rial is lost; therefore, many atoms have unsaturated dan- 
gling bonds. For covalent materials such as Si, the dan- 
gling bond states are located at the center of the band 
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gap, acting as both electron and hole traps, therefore, 
amorphisation of Si is associated with a substantial in- 
crease in electrical resistance, relative to the crystalline 
phaseJH]. For ionic or polar oxides, the dangling bond 
states are very close to the band edges; therefore, amor- 
phisation of oxides usually does not create deep trap lev- 
els, which leads to the remarkable fact that a-TCOs can 
exhibit comparable transport properties relative to their 
crystalline counterparts, an effect that has only recently 
been understood through materials modelling[18, 43-45J. 

The physical and chemical properties of TCO materials 
have been the subject of many excellent reviews. In par- 
ticular, the pioneering experimental work of Hosono et 
al. has been recently reviewed [46], while a more compre- 
hensive review of the field can be found in a themed issue 
of the Material Research Bulletin [47 . The chemistry of 
the ZnO-In 2 03-Sn02 system was the focus of a detailed 
review by Hoel et al.^\. From a theoretical perspective, 
comprehensive overviews focused on ZnO can be found in 
the works of Catlow et aL[49], Lany and Zunger[50] and 
Janotti and Van de Walle[5Tj. In addition, Medvedeva 
and Hettiarachchi recently provided an in-depth analysis 
of the properties of complex TCOs based on trends in 
their physical properties [18]. 

In this Topical Review, we focus on recent advances in 
the theoretical understanding of multi-component TCOs. 
After briefly addressing the binary metal oxide compo- 
nents of the complex TCOs, we address a number of 
highlights from our recent investigations, including: (i) 
the fundamental principles governing the structures of 
multi-component crystalline and amorphous oxide sys- 
tems; (ii) the relationship between composition and opto- 
electronic behaviour; (iii) the origin of high-performance 
for amorphous oxide semiconductors. From these ad- 
vances, future directions in the field of transparent con- 
ducting semiconductors are discussed. 

COMPUTATIONAL METHODS 

The workhorse technique in applied computational ma- 
terial science is Density Functional Theory (DFT)[52"1 
[55] . for which Walter Kohn and John Pople shared the 
Nobel prize in 1998. In contrast to quantum chemi- 
cal methods, which are concerned with the 37V dimen- 
sional many-body wave-function, in DFT it is the elec- 
tron density that is the key quantity to be optimised 
and from which all ground-state properties can be de- 
termined. The result is an ab initio method that can 
be applied to realistic material simulations in solid-state 
physics, organic chemistry, soft matter and related fields. 

Approaches based on DFT are widely available in a 
number of academic and commercial packages, e.g. the 
codes VASP]5Q, CASTEP\§5\, Quantum-Espresso^ 
and FHI-AIMS[57 . The principal difference between 
all implementations is in the choice of the basis set, 



which is used to represent the single-particle Kohn-Sham 
wave-functions, e.g. plane-waves, augmented plane-waves 
with atomic orbitals, gaussian functions, Slater orbitals, 
numerical orbitals. The technical details of these ap- 
proaches have been reviewed in detail elsewhere |58l 159] . 

Evaluation of the total energy, within the DFT for- 
malism, gives access to a range of thermodynamic prop- 
erties including heats of formation, ionisation energies 
and phase stability. Furthermore, forces and stress can 
be calculated routinely, which provides the flexibility for 
a wide range of studies on structure determination, in- 
cluding, for example, global optimisation 60 - 63] ■ The so- 
lution of the Kohn-Sham equations yield the eigenvectors 
and eigenvalues for the respective system, and hence, it 
can be used to understand the underlying bonding me- 
chanics and electrical transport properties. 

The limitations of DFT are concentrated in the 
exchange-correlation functional, which contains the de- 
scription of all quantum mechanical effects for the sys- 
tem of interest. The exact form of the functional remains 
unknown and must be approximated, which is generally 
in the form of local or semi- local functionals, i.e. the lo- 
cal density approximation [51] (LDA) or generalised gra- 
dient approximation 65 (GGA), respectively. Both ap- 
proaches gives rise to an accurate description of the struc- 
tural, electronic and thermodynamic parameters for most 
systems; however, they fail to give a proper description 
of particular properties and systems, in particular, highly 
correlated (strongly interacting) electronic states [55]. 

It has been well demonstrated that both the LDA and 
GGA exchange-correlation functionals underestimate the 
electronic band gaps of insulators and semiconductors 
due to a discontinuity in the wavefunction character at 
the band edges [ST] [55]. For quantitative band gap pre- 
dictions, the use of many-body approaches such as GW 
theory [69] are necessary; however, the use of non-local 
hybrid density functionals offers a means to reduce the 
error within the single-particle formalism. In the hy- 
brid approach, a percentage of exact non-local Hartree- 
Fock exchange is added to a local or semi-local exchange- 
correlation functional. One caveat of this approach is 
that the amount of exact-exchange and spatial screen- 
ing remains material and property dependent [70 . One 
hybrid functional, which has had particular success in 
its application to oxides, is the Heyd-Scuseria-Ernzerhof 
(HSE06) formulation[7T] 172] . in which 25% of exact ex- 
change replaces the short-range GGA exchange potential. 

BINARY OXIDES 

Here, we briefly discuss the bulk properties of the main 
binary components of the multi-ternary systems of inter- 
est, i.e. ZnO, In 2 03, Ga 2 03, and AI2O3. A number 
of the material properties are summarised in Table [TJ in- 
cluding the band gaps, heat of formation, cation effective 
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coordination numbers, and the partial charges resulting 
from Bader analysis |73l I74j of the equilibrium electron 
density. Absolute values of the partial charges derived 
from delocalised electron densities are not unique, and 
are generally smaller than the formal ionic charge [751 176j . 
but relative changes in these values can still be instruc- 
tive. 

ZnO 

ZnO is a material that has a rich solid-state physics and 
chemistry, with a range of applications that exploit its 
electronic, optical and piezoelectric properties [77J. Re- 
cently, magnetic doping of ZnO has been the subject 
of intense study, e.g. Co doped ZnO, where the ferro- 
magnetic coupling between Co d states can be electron 
mediated [78 80 . In its ground-state, ZnO adopts the 
hexagonal wurtzite (wz) crystal structure (space group 
P63ITIC) with lattice constants of Oq XP — 3.250 A and 
Cq XP = 5.207 A, [77] in which the Zn and O ions form 
roughly tetrahedral coordinated environments. The ex- 
perimental electronic band gap of ZnO is 3.44 eV[8T]. 
which is direct in nature. As mentioned in the Introduc- 
tion, the material properties of ZnO have been subject 
of recent review both from an experimental [77j and com- 
putational perspectivegMSH E2 [S3] . 

In 2 3 

The high n-type conductivity of IXL2O3 is exploited in 
many TCO device applications, which has greatly con- 
tributed to the rising cost of In metal [54]. Indium ox- 
ide adopts the body-centred cubic (bcc) bixbyite lattice 
(space group Ia3) with a^ p — 10.117 A[85 , which is 
a defective 2x2x2 superstructure of the fluorite min- 
eral structure. The In cations form roughly octahedral 
structures, which implies an effective coordination num- 
ber close to 6. As mentioned above, R12O3 is widely used 
as a TCO; however, the long-standing band gap problem 
of In203 was only recently resolved from a joint exper- 
imental and theoretical effort [H]. The direct electronic 
band gap of the bulk material is of the order of 2.9 eV[TT]. 
and the large disparity between the electronic and optical 
band gaps (E° pt ^3.6 eV) arises from a combination of 
dipole forbidden optical transitions and conduction band 
occupation [9[ \U \ IT2 l 186] . 

Ga 2 3 

Ga2 03 is a wide band gap metal oxide (E g > 4.5 eV), 
which is used in lasers, as well as a dielectric coating 
in electrical devices 87J. Of the five crystal polymorphs, 
the stable /3 phase (space group C2/m) consists of close 



packed oxide ions, with Ga occupying a combination 
of the tetrahedral and octahedral holes (see the coor- 
dination numbers in Table [I]). The equilibrium lattice 
parameters are a^ p — 12.232 A, 6q XP = 3.041 A and 
c c xp = 5.801 A, with p = 103.7° 88J. Both types of Ga 
ions have a formal oxidation state of III, and their lat- 
tice sites contain small lattice distortions away from the 
ideal polyhedra. A recent hybrid-DFT study has pre- 
dicted that formation of an oxygen vacancy in P-G&2O3 
results in a deep donor centre [89"] . 

A1 2 3 

a-AkiOs is an important dielectric material[90 , with 
a band gap of 9.25 eV[HI], which adopts the hexago- 
nal corundum mineral structure (space group R3/c) with 
a c xp = 4.759 A and c xp = 12.990 A[92], and octahedral 
cation coordination. It is widely used as an insulating 
substrate for thin-film growth. Of particular interest the- 
oretically has been its defect chemistry 93-95] . Recent 
work has focused on the magnetic and spin configuration 
associated with the oxygen interstitial [96 . 

TABLE I. Experimental and calculated material properties 
for a number of binary metal oxides, including the electronic 
band gap, the standard formation enthalpy [97], the cation 
effective coordination number, and the Bader partial charge 
(in electrons). 

Compound E 9 (eV) AHf (eV) Coordination Partial Charge 
"ZnO 3.44|81j -3.63 40 L32 

ln 2 3 2.90 111 -9.60 5.9 - 6.0 2.02 - 2.03 

Ga 2 O a 4.90[98] -11.29 4.0 - 5.8 1.80 - 1.91 
A1 2 3 9.25:;9T| -17.37 6.0 2.62 



RULES GOVERNING STRUCTURE 
FORMATION 

Crystalline multi-component oxides 

The multi-compounds (i?M0 3 ) m (ZnO)„ with R, M = 
In, Ga, and Al (n, m = integers) can be synthesized, for 
example, by solid state chemical reactions using stoichio- 
metric proportions of the ZnO, In 2 03, Ga 2 03, AI2O3 
binary compounds and heating at high temperatures 
(about 1000 — 2000 K). The exact temperature required 
depends on the composition and the desired products [9"9T- 
I103j . A large number of syntheses have been performed; 
and crystal structures were first reported for R = In and 
Ga, with different M. The lowest energy structures iden- 
tified for a number of these compounds are shown in Fig- 
ure [T] 

Except for the spinel structured compound ZnAl204, 
there have been limited reports for the synthesis of sin- 
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gle crystal (Al 2 03)(ZnO)„ compounds for n > 1. The 
Zn-rich alloy was recently grown by pulsed laser deposi- 
tion in the work of Yoshioka et aL |104j . where a super- 
lattice structure was identified, which is consistent with 
the homologous crystal structure models developed for 
the In and Ga compounds, as later verified by electronic 
structure calculations [105] . There is growing evidence 
for the presence of the homologous phase in typical Al- 
doped ZnO samples [1061 1107] . which can contribute to 
lower than expected electron carrier concentrations: ex- 
cess electrons are compensated by additional oxygen in- 
corporation in to the lattice. 



(ln 2 3 ) 1 (ZnO) 5 

'Planar' InO 'Modulated' InO 




FIG. 1. Illustrations of the crystal structures for the planar 
and modulated models of (In203)(ZnO)s and the lowest en- 
ergy configuration of (Ga203)(ZnO)e- Indium polyhedra are 
shaded pink, and Ga polyhedra are shaded green, while oxy- 
gen and zinc atoms are coloured red and grey, respectively. 
Visualisation was performed with the VESTA package [108] . 



(IriiOz) m (ZnO) n and (Ga 2 03) m (ZnO)„ 

IZO was first synthesized by Kasper[109] for m = 
l,rt = 2 — 7, who showed that the compounds adopt 
layered structures and crystallise in rhombohedral or 
hexagonal crystal lattices. Two decades later, Can- 
nard and Tilley confirmed Kasper's results using X- 



ray diffraction (XRD) and high-resolution electron mi- 
croscopy (HREM) exp eriments |110j . Recently, Kimizuka 
et a/.[99l 1111] proposed that IZO is isostructural with 
LuFe0 3 (ZnO)„ for m = 1 [ITS]. The crystal has the 
space group R3m (rhombohedral lattice) when n is odd 
and PQ^/mmc (hexagonal lattice) when n is even. The 
In atoms form an In0 2 octahedron layer, which are in- 
terconnected by (InZn„)0„+i layers. The In atoms in 
the (InZn n )O n+ i unit were assumed to be randomly 
distributed on the metal sites, forming 5-fold trigonal 
bipyramidal polvhedra |112] , while the Zn atoms were as- 
sumed to be located in tetrahedral sites surrounded by 
four O atoms. The formation of an In02 octahedron layer 
was also confirmed by atomic-resolution Z-contrast ex- 
periments for In 2 03(ZnO) ra [TT3]. In contrast to the XRD 
studies,[SIS[OT[m! recent HREM experiments [TTMIIH] 
have indicated that the In atoms in the (InZn„)0„+i 
layers form ordered modulated structures with a zig-zag 
shape. 

So far, the synthesis of GZO has been reported only 
for m = n = 1, PZH QH] m = l,n = 6,[M] m = 
l,n = 7,8,9,16,[M] and to = l,n = 9,13. QUI] To 
our knowledge, there is no report for m > 1. Ex- 
perimental characterisation using XRD has found that 
GZO crystallizes in orthorhombic structures with space 
group Cmcm for n > 1, and the 'normal' cubic spinel 
structure for m — n — 1. The internal structural pa- 
rameters (atomic positions) have been reported only for 
m = l,n = 1,6 |103j : however, it is important to note 
that the atomic number of Ga and Zn differs only by one 
unit, and hence, XRD alone cannot clearly distinguish 
between the Ga and Zn site occupation, i.e. it was as- 
sumed that Ga atoms occupy specific cation sites, while 
Zn occupy the remaining sites. The pioneering work re- 
ported by Michiue et al.{l03 provided a basis for the un- 
derstanding of the atomic structure of GZO compounds; 
however, a number of questions remained to be solved at 
that time. 

Recently, Da Silva et aL |119l 1120) employed first- 
principles computational tools to address the mechanisms 
that lead to the formation of the layered IZO structures 
with the In-modulations in the (InZn„)O n +i layers and 
the atomic distribution of the Ga and Zn cations in the 
GZO lattice, which forms natural grain-boundaries. Be- 
low, we will summarize the most important results, which 
can be used as guidelines for a better understanding of 
these materials and other related oxide compounds. 

The preference of In to form octahedron networks: 
The In and O atoms in the bixbyite In 2 03 structure 
are 6- and 4- fold coordinated. |121| respectively, whereas 
both Zn and O atoms are 4-fold coordinated in the wz 
structure. In contrast to the In atoms, the Ga atoms 
in the monoclinic /3-Ga 2 03 structure are separated into 
two groups composed of four distorted octahedron and 
four tetrahedron, while there are eight 3-fold and four 
4-fold O atoms. In particular, a pair of octahedrons are 
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surrounded by 6 tetrahedrons, which can be seen as a ba- 
sic motif in /3-Ga2 03. For the simplest multi-compound 
structure, i.e. the high symmetry i? 2 Zn04 (R = In, Ga) 
spinel structure may also exist, in which the In, Zn, and 

atoms are 6-, 4-, and 4-fold coordinated, respectively. 
The coordination environments observed for the In, Ga, 
Zn, and O atoms are related to the fact that In, Ga, and 
Zn have three, three, and two valence electrons, respec- 
tively, whereas oxygen always has a valence of —2, so the 
coordination environments are determined such that the 
octet rule for local charge neutrality is satisfied. Thus, 
it can be expected that In should prefer to form octa- 
hedrons, whereas Ga can form both octahedrons as well 
as tetrahedrons, even though both atoms have the same 
formal valence. To test this hypothesis, we investigate 
several IZO and GZO model structures. 

Two types of 6-fold coordinated InC>2 layers were calcu- 
lated in IZO, namely, a non-octahedron structure and an 
octahedron structure. The octahedron structure is about 

1 eV/In02 lower in energy than the non-octahedron 
structure. This is consistent with the fact that the oc- 
tahedron structure maximizes the atomic separation be- 
tween the negatively charged O atoms. In the octehedron 
structure, the In— O distances in the In02 layers are in 
the range of 2.20 - 2.30 A, while the angle of O-In-0 
deviates by 1 — 5° from the ideal value of 180°. The 
smallest deviations (1°) occur when O atoms in the cor- 
ner of the octahedron are shared with Zn atoms. Thus, 
the binding of Zn to O at the corner of the L1O2 octahe- 
dron layers stabilises the formation of an almost perfect 
In-0 octahedron, whereas binding to In induces a larger 
distortion. These findings are consistent with the for- 
mation of distorted and perfect octahedron structures in 
the bixbyite In203 and the spinel In2Zn04 structures, re- 
spectively. The observation that the In and O atoms in 
the In02 layer form an octahedron structure is consistent 
with experimental observations. |991 1111H113 1 1122] 

In contrast to the In atoms in IZO which is 6-fold co- 
ordinated, Ga atoms may exist at low coordination sites. 
For example, in the GZO structure proposed by Michiue 
et al. |103| there are two cation sites (called M8 and Ml 
in Ref. ll03p . which are 5- and 4-fold coordinated and were 
assumed to partially occupied by Zn atoms. We showed 
that the M8 sites are empty and the Ml sites are oc- 
cupied by Ga atoms, while the remaining Ga atoms are 
distributed in 6- fold like sites. Similar to the fi-Q^Oz 
structure, both Ga sites satisfy the octet rule. Thus, 
these results show clearly that In atoms have a stronger 
preference for octahedron sites in IZO than Ga atoms in 
GZO, which explains the formation of planar In02 layers 
in IZO, whereas it is absent in GZO. 

Inversion domain boundary formation in (RZn n )O n +\ 
layer: As discussed above, the In atoms have a strong 
preference for octahedron sites, while the Ga atoms can 
occupy 6- and 4-fold sites. Thus, once the 6-fold In sites 
in the In02 layer in IZO and the 4-fold Ga sites in GZO 



are occupied, the remaining In and Ga atoms must be 
distributed in the wz-like ZnO structure; however, all 
cation sites in the wurtzite structure are 4-fold and not 
6-fold. Thus, changes must occur in the (i?Zn„)0„ + i 
layers in order to accommodate both In and Ga atoms 
and satisfy the electronic octet rule. 

Furthermore, a number of structural features of IZO 
and GZO also contribute to the formation of an in- 
version domain boundary (IDB). For example, in IZO, 
the In02 octahedron layers are connected by the wz-like 
(InZn„)0„+i layers, in which the O atoms at the corner 
of the octahedrons are connected to three In atoms in 
In02 and one atom in the (InZn„)0„+i layer. Thus, 
the polarity at the bottom and top of the repeating 
unit is reversed, and hence, an IDB must exist in the 
(InZn„)0„+i layer to reverse this polarity. An In02-like 
layer is not present in GZO; however, the distribution of 
the Ga atoms requires 6-fold like sites, which can be ob- 
tained through the formation of IDB in GZO. For both 
IZO and GZO cases, the IDB is located on the In and Ga 
atoms, which form 5-fold trigonal bipyramid structures 
with the surrounding O atoms in the (-RZn„)O n +i layers. 

Stacking fault in the (RZn n )O n +\ layer: It is impor- 
tant to realize that the formation of the IDB does not 
occur alone, i.e. 5- fold trigonal bipyramid sites cannot 
be created without disrupting the stacking sequence of 
the wz-like lattice. In order to preserve the hexagonality 
of wurtzite and without destroying the stable In02 octa- 
hedron layer and the formation of the 4-fold Ga sites, one 
or more stacking faults must exist in the (i?Zn„)0„+i 
layers. This has important consequences for determin- 
ing the stable configurations of IZO. For example, the 
conventional hexagonal unit cells are composed of two, 
or three, In02 layers separated by an equal number of 
(InZn„)0„ + i layers in which n is even or odd. This leads 
to a hexagonal primitive unit cell for n even and rhom- 
bohedral (or monoclinic) for n odd. 

Modulation of R atoms in the (RZn n )O n+ i layer: As 
previously mentioned, to obey the material stoichiomc- 
try, a number of In and Ga atoms must adopt trigonal 
bipyramid structures with the surrounding O atoms; at 
low temperature, these motifs should form ordered con- 
figurations. The ionic radii of the In, Ga, and Zn atoms, 
i.e. 0.80 A for In(III), 0.55 A for Ga(III), and 0.60 A for 
Zn(lI). [T2"5] play a decisive role in determining the or- 
dered structure that minimizes the internal strain of the 
IZO and GZO structures. 

For example, model calculations demonstrated that 
the in-plane lattice constant of InM03 in the hexago- 
nal InFeOs structure |122j is 8.4% larger than the in- 
plane lattice constant of ZnO in the wz structure. Thus, 
the In atoms must assume an ordered distribution in the 
(InZn„)O r ,,+i layers to minimize the strain energy. Cal- 
culations for several different atomic configurations have 
identified that the In atoms are distributed in the cation 
sites following a zig-zag modulated structure in IZO, in 
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which the modulation period is proportional to the num- 
ber of ZnO units. Similarly, the Ga atoms are distributed 
in GZO to form a modulation, which match the one ob- 
served in IZO. 

The modulation is composed of In/Ga and Zn atoms, 
which form trigonal bipyramidal structures. These re- 
sults are consistent with experimental HREM studies, 
[11514118] which observed clearly the formation of In- 
modulated structures in the (InZn n ) n+ i layers for n > 5 
ZnO units and m = l |116j . 

Electronic octet rule for the O in the (RZn n )O n+ i 
layers: The minimum energy structure for i?Zn n O n +i 
obeys the electronic octet rule (local charge neutrality), 
in which the presence of 5-fold Zn sites is fundamental 
for both IZO and GZO. The following environments were 
identified in IZO as fundamental for the electronic octet 
rule: the O atoms are surrounded by 4xZn(4) atoms, or 
2xZn(4) + lxln(5) + lxZn(5), or 2xln(5) + 2xZn(5), 
where the numbers in parentheses indicate the integral 
coordination of the In and Zn atoms, so that each O atom 
always receives two electrons from neighbouring cations. 
Similar coordination environments are observed in GZO 
beyond of the 4-fold Ga sites. Furthermore, the octet 
rule is also satisfied for oxygen pairs along the in-plane 
In and Zn rows, e.g. 2 x Zn(4) + 3 x Zn(5) + 3 x In(5), 
and for O atoms in the In02 layers. Therefore, the local 
electrostatic environment plays a fundamental role in the 
stability of the (i?203) m (ZnO)„ compounds. 

(InM 111 '0 3 )m(ZnO)„ 

One important requirement for technological applica- 
tions of IZO is complete control of the electronic proper- 
ties, e.g. band gap, and hence, several studies have been 
performed to obtain a procedure to tune the band gap 
of IZO for specific applications. One possible approach 
is the replacement of a percentage of In atoms by other 
trivalent species, e.g. Ga or Al atoms. Therefore, it is 
very important to understand the consequences of such 
substitutions on the IZO atomic structure. We observed 
that the ionic size of the In(III) and Ga(III) cations rel- 
ative to Zn(II) plays an important role in the formation 
of the M-modulations. 

In our calculations, we assumed that half of the 
In atoms were replaced by Ga or Al atoms to form 
(InGa03) m (ZnO)„ and (InA103) m (ZnO)„, respectively. 
We found that the In atoms always form the In02 layer, 
which is expected based on the high stability of In octahe- 
dron structures. Replacement of the In atoms in the In02 
layer by Ga or Al atoms and put In atoms in the ZnO 
layer, is highly energetically unfavourable. Thus, Ga and 
Al should replace the modulated In atoms in the ZnO lay- 
ers. We also demonstrated that the in-plane lattice con- 
stant of InA/03 in the hexagonal InFe03 structure |122j 
is larger by 8.4%, 2.4%, and 0.0% for M = In, Ga, and 



Al, respectively, compared with the lattice constant of 
wz-ZnO. Therefore, for M — In, a large strain exists in 
the ZnO and InO layers if the In atoms form a planar 
layer. The in-plane strain decreases for M = Ga and is 
almost zero for M = Al. Therefore, it is probable that 
the Al atoms are randomly distributed in IAZO, while in 
IGZO, modulated Ga ions might still be present. 

Amorphous multi-component oxides 

In contrast to the high temperature associated with the 
synthesis of crystalline multi-component oxides, amor- 
phous thin films can be grown by low temperature pro- 
cessing. Following an initial investigation by Hosono et 
al. on amorphous ternary oxides [T2"4"] . the breakthrough 
for a-TCOs came from the report by Nomura et aL[3"5] 
that the electron mobility of amorphous IGZO thin films 
could exceed 10 cm 2 V _1 s _1 , which is an order of magni- 
tude greater than amorphous Si. IGZO has since become 
the most widely studied amorphous TCO, and has found 
application as an electron-injection-layer for solid-state 
lighting devices [31 [HS] EH] . 

To simulate these low symmetry structures, a pseudo- 
cubic supercell of the crystalline structures was gener- 
ated, which was subject to high temperatures under ah 
initio molecular dynamics until a solid to liquid phase 
transition was observed. Then, the liquid oxide was 
quenched, and the resulting structures were subject to 
local optimisation of both the internal positions and lat- 
tice constant. The full details of the approach employed 
can be found in Reference H3] A snapshot of an amor- 
phous structure is shown in Figure [2j 

As previously shown, the crystal structures of the In 
and Zn containing oxides are defined by mixed In-0 and 
Zn-0 polyhedra (bipyramids and tetrahedra) present be- 
tween In02 octahedral planes. In the amorphous phases 
of IZO, IGZO and IAZO, the same local coordination 
motifs are maintained. To compare the local order in 
both phases, we have calculated both the effective coor- 
dination number [127] 1128] and effective charge distribu- 
tion for each ion[43]. Even in the amorphous phase, the 
strong charge transfer between the metal cations and oxy- 
gen is preserved, driving the similar local coordination. 
While this is to be expected due to the large electro- 
negativity of oxygen, it is in contrast to covalent semi- 
conductors, where amorphisation can lead to a large and 
significant local bond rearrangements that perturb the 
nearest neighbour environments [129J. 

The total energy difference between the crystalline and 
amorphous TCOs is of the order of 200 meV per f.u. in 
each system studied, which is expected based on the in- 
creased strain due to the disordered packing of the cation 
centred polyhedra on amorphisation. Despite their high 
energy, the a-TCOs remain metastable experimentally at 
standard temperatures as they do not have sufficient en- 
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ergy to overcome the large kinetic barriers required for 
recrystallization. 




(lnAI0 3 )(ZnO) 1 



FIG. 2. Conduction band charge density isosurfaces 
(coloured turquoise) for two crystalline and one amorphous 
oxide material. The oxygen ions are coloured red. 



OPTICAL AND ELECTRONIC PROPERTIES 

The performance of electroceramics for TCO appli- 
cations can be quantified by their electrical conductiv- 
ity and optical absorption spectra. The former is re- 
lated phenomenologically to the position of the elec- 
tronic states on an absolute energy scale (band align- 
ment), which influences the material's doping limits, 
while the latter depends on the magnitude of the funda- 
mental band gap, dipole selection rules for optical transi- 
tions, and carrier concentrations (e.g. the Moss-Burstein 
effect pPlH5ip . 

Beyond the bulk material properties, the role of de- 
fects and non-stoichiometry in these oxides is also im- 
portant in determining their performance for TCO ap- 
plications. The theoretical framework for the accu- 
rate modelling of point defect formation exists, and 
the resulting free energies of formation for each defect 
charge state can be solved self-consistently to produce 
temperature dependant equilibrium defect and carrier 
concentrations [T^l I133J . Quantitative defect modelling 
techniques have not yet been applied to the case of the 
complex multi-component systems. The difficulty lies in 



the large number of inequivalent lattice positions, which 
must be calculated independently, as well as the lower lat- 
tice symmetry in comparison to the binary components; 
however, with the increase of high performance comput- 
ing resources to the peta-scale, these issues will soon be 
overcome. Here, we can rationalise the material perfor- 
mance based on the bulk electronic changes observed in 
the multi-component systems. 

Natural band offsets 

The band offsets between semiconductors is one of the 
most fundamental properties for device physics [I34j . The 
absolute band edge positions of a material have a role in 
determining intrinsic doping limits. I136j as well as 
their electrochemical properties [137] and interfacial elec- 
tron transport. Experimentally, X-ray Photoemission 
Spectroscopy (XPS) can be used to measure the posi- 
tion of atomic-like core states, which provide a reference 
for the alignment between the valence bands of two or 
more systems [138| . A theoretical approach has been de- 
veloped in the same spirit, from which the 'natural' va- 
lence band offset of two compounds at their respective 
equilibrium positions can be calculated |139| 1140] ; this 
approach has recently been revised to take into account 
the deformation of the core states with respect to vol- 
ume changes [141j . The conduction band offsets can be 
inferred by adding the fundamental band gap to the com- 
puted valence band offset. 

The natural band offsets derived for a number of oxide 
materials of interest are shown in Figure [3j The results 
match the recent alignment of the ultraviolet and inverse 
photoemission spectra for ZnO and In9 0^ [T4"2] . The ob- 
served trends are determined by a combination of shallow 
cation d states (e.g. Cd 4d states raise the valence band 
of CdO), the degree of ionicity (e.g. efficient electron 
transfer produces a high electrostatic potential that low- 
ers the valence band of MgO), and crystal structure (e.g. 
low coordination oxygen sites raise the valence band for 
amorphous IZO). 

It can be observed from Figure [3] that I112O3 has the 
lowest conduction band level (highest electron affinity) 
of all the transparent oxides studied, which correlates 
with its ability to sustain high electron carrier concen- 
trations that result in degenerate electron conduction. 
Indeed, recent experimental work has shown that the 
charge neutrality level, i.e. the Fermi energy at which 
oppositely charged defects can form spontaneously, of 
In2 03 lies well above the conduction band[S], which has 
been supported by calculations at the level of many-body 
perturbation theory 143 ; further detail and implications 
of the charge neutrality level are discussed in work by 
King in this themed issue |144] . While ZnO is a good 
n-type semiconductor, the higher conduction band level 
limits its ability to be heavily electron doped. For large 
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gap oxides such as AI2O3 and Ga 2 03, electron doping 
will be even more difficult, with high degrees of ionic 
charge compensation expected to occur [145]. These ef- 
fects have been summarised in a series of doping limit 
rules for semiconductors [T^fl 11361 fl46| . 

Forbidden optical transitions and carrier localisation 
in (In 2 03)(ZnO) n 

The principal reason why In 2 03 remains a good TCO, 
despite its relatively small band gap, is that the band 
edge optical transitions are symmetry forbidden: it 
combines a small fundamental gap with a large op- 
tical band gap. Can this also be the case for the 
ternary In-containing compounds? For the homologous 
(In 2 03)(ZnO)„ series, we have shown that the lowest en- 
ergy optical transitions from the valence to conduction 
bands are dipolc allowed for all values of n, as the top of 
the valence band state has large contribution from ZnO 
|148j . This explains the anomaly of why the apparent 
optical gap of ln 2 3 (ca. 3.6 eV)[7J reduces to below 3.2 
eV on the addition of ZnO to form IZO [51 1551157]. 

For the IZO compounds, while the top of the valence 
band is more localised on the Zn-0 network, the conduc- 
tion band is localised on the In-0 network, which follows 
from the band alignment of the binary materials (Fig- 
ure pj). For n = 1, the localisation is not significant as 
there are two mixed cation layers present between the 
In octahedra, and these contain an equal number of In 
and Zn ions. For higher values of n, the number of in- 
termixed layers increases, which raises the Zn:In ratio, 
and the band edge wavefunctions become more localised 
as a result. With the increased Zn concentration, an 
effective superlattice is formed and the charge carriers 
become confined in the In0 2 octahedron network. The 
modulation of In ions through the ZnO layers makes a sig- 
nificant additional contribution to the conduction band, 
producing a modulated electron density distribution, as 
illustrated in Figure [2] The predicted carrier localization 
correlates well with the minimum of resistivity previously 
found for IZO samples with low Zn:In ratio, i.e. for the 
range < n < 1 [5J [55J [57] . 

High performance of amorphous oxides 

For the crystalline TCO compounds, the band gaps 
follow the clear trend of the group 13 binary oxides, in- 
creasing from In to Ga to Alj35J. For the amorphous 
compounds, the same band gap trend is observed, but 
the magnitudes are lower compared to the crystalline 
structures. The differences arise from the series of bands 
situated at the top of the valence band. In the crys- 
talline oxides, the lowest binding energy valence bands 
are dominated by O 2p states, and this is unchanged 



for the amorphous systems; however, the variation in 
electrostatic potential around each of the oxygen sites 
in the amorphous structures results in the spatial local- 
isation of the valence wavefunctions. The contribution 
of the O 2p valence band to the sub band gap states 
has been confirmed through hard X-ray photoemission 
measurements [126], and is consistent with the band gap 
decrease observed for a-TCOsj35j. As these states are 
concentrated in the Zn-rich regions (due to the lower 
binding energy shallow-core Zn 3d states), their contri- 
bution to the absorption of visible light is not significant 
because the conduction band is more localised on the 
group 13 cation side, and the resulting optical transition 
probabilities are low. As a result, the amorphous materi- 
als retain sufficient transparency in the visible region to 
be viable for TCO applications. 

The electronic density of states predicted for the amor- 
phous oxides has important consequences. Crucially for 
n-type conductivity, the conduction band remains highly 
delocalised, which can be understood from the overlap 
between the disperse cation s orbitals and the isotropic 
nature of the cation s - O 2s bonding[551 |4"H 1149] : the 
overlap is largely unaffected by the long range structural 
disorder, which is illustrated in Figure [2] 

As expected from their larger In— O bond length com- 
pared to the other cations, In makes the largest contri- 
bution to the antibonding conduction band state. For 
IZO, two thirds of the cations are In, while for IGZO 
and IAZO, only one third In are present. As stated ear- 
lier, even for the lower concentration materials such as 
IAZO, the In— O polyhedra form a combination of edge 
and corner sharing networks ensuring that the conduc- 
tion band is appropriately spread over the entire amor- 
phous network. This also explains why the amorphous, 
low indium concentration TCOs can exhibit good electri- 
cal transport properties: in the crystalline compounds, 
the In— O networks become confined to two dimensional 
planes separated by the Zn and group 13 cations, while 
in the amorphous compounds, the indium ions become 
evenly distributed and hence offer improved isotropic car- 
rier mobility. 

While a-TCOs exhibit excellent n-type conductivity, 
p-type behaviour is not expected to be present due to 
the strong localisation of the upper valence band. The 
majority of known p-type oxide semiconductors such 
CuA10 2 and SrCu 2 2 are derived from Cu 2 0, where 
electron deficiency (hole formation) is facilitated by Cu 
oxidation, with carrier mobility along the linear O-Cu-0 
chains. Here, electrical carrier (polaron) mobility will 
be inhibited by the structural disorder on amorphisa- 
tion. The absence of p-type response has been reported 
experimentally[126 . The discovery of a high perfor- 
mance p-type TCO material, either amorphous or crys- 
talline, remains an open challenge in this field. 
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FIG. 3. Natural band alignment between In203, ZnO and related oxides, relative to that of wurtzite ZnO, and compiled 
from References I43l IT4TI 1 1471 and 11481 The differences between the electronic and optical band gaps of In2 03 and amorphous 
(ZnO)(In203) are emphasised. The placement of the charge neutrality level is based on an explicit calculation for In203. 



DISCUSSION AND OUTLOOK 

It is clear that the formation of multi-component ox- 
ides offers a way to tune the physicochemical properties 
associated with transparency and conduction. We have 
highlighted a number of rules that determine the struc- 
ture formation of these complex materials. While the 
crystal structures of the binary oxide components are well 
described by the established principles of Pauling 150, 
I151j . for the multi-cation systems, features such as polar- 
ity inversion, stacking faults and partial site occupation 
arise. In terms of the electronic properties, for these n- 
type materials, the key attribute is a delocalised conduc- 
tion state, which can be maintained when a significant 
amount of In is present. The mixing of Ga and Al serves 
to both increase the band gap and the crystal stability, 
which can provide a means to control the stoichiomctry 
and optical transmission, e.g. to provide a blue-shift in 
the optical band gap suitable for UV applications. 

Despite these advances, a number of issues are still out- 
standing in the field of multi-component TCOs. Even in 
the 'simpler' binary metal oxide components, the origin 
and nature of conductivity remains a contentious issue. 
One common feature is that oxygen vacancies appear to 
be deep donor centres, with the neutral defects being 
stabilized by the strong Madelung potential of the lat- 
tice; conversely, metal interstitials generally act as shal- 
low donors, but have higher formation energies. The role 
of hydrogen as an extrinsic donor impurity has been pro- 
posed for a number of TCO materials }5 1 U 1 52l4154j . but it 
is highly unlikely to be the primary electron source. For- 



tunately, atomistic simulation techniques are becoming 
increasingly quantitative, and we can now predict defect 
formation energies with high accuracy, especially when 
temperature dependent contributions to the free energy 
are included. It is therefore anticipated that our under- 
standing of the binary oxide systems will be solved in the 
near future, and that this understanding will transfer to 
the more complex multi-component systems. 

One real limitation to transparent electronics is the 
absence of p-type TCO materials, which are needed to 
form p — n junctions for diodes and transistors. It ap- 
pears that the performance of Cu-based compounds will 
not be sufficient for commercial applications. An alter- 
native route to the use of Cu(I) is to form oxides of 
cations with filled s bands such as Sn(II) or Bi(III) [155J . 
Both SnO and BiV04 have attracted recent attention for 
their p-type conductivitv |1561 1157] . Moving away from 
pure oxide materials, forming solid-solutions with chalco- 
genides (S, Se and Te) may offer a way to increase hole 
stability by raising the valence band energy. In addi- 
tion, the formation of hybrid materials containing n-type 
inorganic networks and p-type organic networks is also 
feasible [l58j . and we have performed exploratory work 
in this area[T59 l [T60 ] . 



CONCLUSION 

We have demonstrated the success of computational 
modelling techniques applied to the complex case of 
multi-component transparent conducting oxides. In the 
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last decade, our understanding of the structural, elec- 
tronic and optical properties of these materials has in- 
creased significantly. A number of fundamental 'rules' 
have been established relating to crystal structure for- 
mation, as well as the conductivity trends, which can 
provide a framework to enable, for example, a high- 
throughput computational screening of candidate com- 
pounds and stoichiometries. We note that despite these 
advances, many open questions remain in the field, espe- 
cially concerning the microscopic origin of conductivity 
and the role of extended defects, such as dislocations and 
grain boundaries, which ensure that the field will remain 
active for the near future. 
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